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Abstract 

We derive a path-integral description of the vortex state of a fermionic superfluid in the crossover 
region between the molecular condensate (BEC) regime and the Cooper pairing (BCS) regime. 
This path-integral formalism, supplemented by a suitable choice for the saddle point value of the 
pairing field in the presence of a vortex, offers a unified description that encompasses both the BEC 
and BCS limits. The vortex core size is studied as a function of the tunable interaction strength 
between the fermionic atoms. We find that in the BEC regime, the core size is determined by 
the molecular healing length, whereas in the BCS regime, the core size is proportional only to the 
Fermi wave length. The observation of such quantized vortices in dilute Fermi gases would provide 
an unambiguous proof of the realization of superfluidity in these gases. 
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I. INTRODUCTION 

The realization of superfluidity in dilute fermionic gases has opened up a new av- 
enue for the investigation of fermionic quantum systems. In such gases, the presence of 
Feshbach resonances in the interatomic scattering allows to tune the interaction strength 
from strongly repulsive to strongly attractive 2\. In the limit of weak repulsive interactions, 
deeply bound molecular bosons are formed that Bose condense [l| (BEC). In the limit of 
weak attractive interactions, a Bardeen-Cooper-Schrieffer (BCS) fermionic superfluid arises 
P). The crossover between both regimes has recently attracted a great deal of theoretical 
and experimental interest. 

One of the hallmarks of superfluidity, be it bosonic or fermionic, is the presence of quan- 
tized vortices. Whereas vortices are well understood both in the BEC and the BCS limits, 
it is not clear how the characteristics of the vortex (such as the core size) behave in the 
crossover regime. Bulgac and Yu have extended density functional theory to superfluid 
fermion systems |4|, and studied vortex states in the BCS regime within their superfluid 
local density approximation (SLDA) (sj). They found that in the BCS regime vortices give 
rise to a depletion in the fermion density. These results confirm earlier calculations based 
on the Bogoliubov-de Gennes theory 

The goal of this communication is to develop a path-integral treatment suited to describe 
vortices in the BCS, BEC, and crossover regimes. With this treatment, we investigate how 
the vortex core size and the fermionic density depletion at the core change when the fermionic 
superfluid is brought from the BEC to the BCS regime. On the BCS side of the Feshbach 
resonance, we compare our results to those obtained with the SLDA treatment Q|. 

A path-integral treatment for the ground state of the fermionic superfluid was developed 
by Sa de Melo, Randeria and Engelbrecht j?l 1^- Their formalism provides a unified de- 
scription of the BEC, BCS, and crossover regimes and predicts a smooth cross-over for the 
critical temperature, pairing gap, and chemical potential. They consider a homogeneous 
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Fermi gas of atoms determined by the action functional 



S 1 = dr dx J2a 



^ ( £ - ^ v x - H ^ 



+ / dT / dx ^V'x,r,T'0x,r,4'0x,r,4'0x,r,T- 



(1) 



In this expression ^x,r,cr and ^x.t.o- are the Grassmann variables describing the fermionic 
degrees of freedom, where x is the position vector, r the imaginary time and a =|, | denotes 
the two hyperfine spin states present in the Fermi gas. The chemical potential is denoted 
by /i, and (3 = 1/(/cbT) is the inverse temperature. The interaction between the fermionic 
atoms only takes place between atoms in different hyperfine spin states, and is described 
by a contact interaction characterized by the renormalized strength g. The partition sum is 
given by the functional integral over the Grassmann variables, 



Z = J V[lpx,T,a, 1px,r,a] exp {Si} 



(2) 



To unravel the product of four Grassmann variables, the Hubbard-Stratonovic transforma- 
tion is performed. This transformation introduces the bosonic Hubbard-Stratonovic fields 
A XjT and A XjT such that 



Z= £>[A XjT , A XjT ] / V[ip x>Tia ,^ X:Tja \exp{-S 2 } . 



(3) 



with 



p 
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A x , T A XjT 
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(4) 



The action functional S 2 is quadratic in the Grassmann variables so that the functional 
integration over these variables can in principle be evaluated. 



II. SADDLE POINT FOR THE VORTEX STATE 



In practice, to perform the functional integration over Grassmann variables, one has to 
choose a saddle point value for the fields A XjT , A XiT . To describe the ground state, Sa de 
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Melo et al. [7] suitably choose a uniform constant saddle point (A X)T )* = A X)T = |A|. After 
doing this, the Grassmann variables can be integrated out straightforwardly, resulting in an 
effective saddle-point action 

S eS = Jdr J dx I ~ti [In (-GT 1 )] - J-L L , (5) 

where G _1 is the inverse Nambu propagator, for the ground state: 



Aground state 
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or 



0"3 



2m 



(6) 



and the <jj are Pauli matrices. The saddle point equation SS e g/8 |A| = then leads to the 
familiar gap equation for the fermionic superfluid The chemical potential \x is fixed by 
the fermion density. At finite temperatures, fluctuations around the saddle point value can 
be taken into account to improve the theory and find the critical temperature [3, [| • 
To investigate the vortex state, we propose to use a different saddle point, and set 



(A X)T )* = A X;T = |A r | exp(i6 



(7) 



where 6 is the angle around the vortex line and r is the distance to the vortex line. That 
is, r and 9 are the radial and angular coordinates if one chooses cylindrical coordinates 
x = (r, 6, z) such that the z-axis lies along the (straight) vortex line. This particular choice 
of the saddle point value lies at the core of the present treatment. With this choice for 
the saddle point value, the integration over Grassmann variables leads again to an effective 
action of the form (jHJ), but with a different result for G _1 . We find: 



^Vortex 



= CTO 
+ 0"3 



d i 
dr 2mr 
1 



a 1 1 A r 



^V 2 

2m x 



(8) 



The density of paired atoms near the vortex core (|A r | as a function of r) can in principle be 
derived from the saddle-point equation 5S e s/S \A r \ =0 (although in practice the presence of 
the spatial derivatives inhibits straightforward calculation). However, the result no longer 
contains explicit information on the density of fermionic atoms. The results of Bulgac et al. 
interestingly show that whereas |A r | tends to zero at the vortex core, the fermionic density 
needs not go to zero. To study the density of fermionic atoms, expression (JHJ) is not suited. 
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III. EXPLICATING THE FERMIONIC DENSITY 



In order to introduce the fermionic atom density in the path-integral expressions, a trans- 
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formation was proposed by De Palo et al. |l(| based on the identity 



p 

C = J V \p™, p x , r ] exp -Jdrjdx ip™ 

L o 

X (Px,r - ^x,t,T^x,t,T ~ ^x,T,iV>x,T,i)] 



(9) 
(10) 



where C is a constant c-number and p^,p X)T are bosonic fields. This generalized delta- 
function expression identifies p XjT with $^ .V'x 1 t,<tV'x ) t,<t> the fermionic density. 

Multiplying the partition sum (j3J) with the constant f(10j) . and approximating the 
Hubbard-Stratonovich fields A XjT , A X;T by the vortex saddle point (J7J) leads to the following 
saddle-point approximation for the partition sum 



Z sp = J V [ip x , T ,a,ip x ,T,a] Jv[p^,Px,t] exp{-S sp } 



(11) 



with 
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IAJ 



(12) 



To get rid of the path integration over p^ff , p XiT in ([lip, we again use a saddle-point approach 
and set these fields equal to p^f = p^ s , p x>T = p r . This implies that the fermion density 
only depends on the distance r from the vortex core (and not on the 6 or z coordinate). 
From (|T2|) it is clear that the auxiliary field ip^ s is related to the chemical potential. It is 
useful to introduce [ll3| 



Cr 



ip" S + p 



ymr 1 



(13) 



On the level of the fermionic degrees of freedom, ( r acts as a one-body potential combining 
the auxiliary field, the chemical potential and the angular momentum barrier. Alternatively, 
it can be interpreted as a local chemical potential. 



To perform the integration over the Grassmann variables we assume that ( r and the 
pairing field A r vary slowly in comparison with the relevant fermion frequencies. After 
integration over ■0 x ,t,<tj V'x.r.o-j we then obtain the result Z sp oc exp {S' eS } with 

12 ' h 2 



I At, 
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(r-fl 



■smr' 



Pr 



+2 



2cosh(l V /(£- C ,) 2 + |A r | 2 j]}. 



(14) 



IV. PROPERTIES OF THE VORTEX CORE 



The first saddle point equation is 



9S' eS 



dp r 



<^> Cr = P - 



(15) 



This locks the local chemical potential ( r as a function of the overall chemical potential and 
the energy barrier of the vortex flow, as one would expect in the Thomas-Fermi approxima- 
tion in the molecular limit. The parameter [i is determined by fixing the total number of 
particles (or the density far away from the vortex core). Differentiating S' eS with respect to 
|A r | we find the familiar gap equation, but with a local chemical potential determined by 

Cr'- 
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with Ek 



- I dkk 2 



Aanh(/3£ fc /2) 
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(16) 



(P - C,) 2 + A?. 



Finally, the third saddle point equation dS' eS /d( r 
density near the vortex core 



allows to calculate the fermionic 



Pr 



| / dk k 2 



k 2 £ r 
El 



tanh(/3E fc /2) 



(17) 



In ([16 )1 . ([17 )1 . units are such that wave numbers are expressed in Fermi wave numbers and 
energies in Fermi energies. 

We investigate the solutions of these equations in the zero temperature limit (j3 — ► oo), 
so that we may assume that fluctuation corrections around the saddle point values are not 
important. The result for the pair density and the fermion density are shown in figure 1. 
Both the pair density and the fermion density go to zero in a region with spacial dimensions 




FIG. 1: The pair density relative to the bulk pair density (left axis, dashed curves) and the 
fermionic atom density relative to its bulk value (right axis, full curves) are shown as a function of 
the distance to the vortex line, for two different values of the interaction strength. The depletion 
of the fermion density at the core is overestimated by the current approach, as discussed in the 
text. The thin part of the curves indicate the region where the current approach is unreliable. 



of the order of l/kp around the vortex line (r = 0). In the BEC limit, the pair density and 
the fermionic atom density become zero at the same distance from the vortex line, so that 
no atoms are seen inside the core. This corresponds to the expectation for a molecular BEC. 
In the BCS limit, fermionic atoms can penetrate into the region where the pair density is 
zero. These results qualitatively agree with those of Bulgac et al. in the BCS region, but 
quantitatively they are quite different: in refs. 0,|f| only a small depression in the fermionic 
density is found on the BCS side of the Feshbach resonance. What can be the reason for the 
discrepancy ? Most likely the assumption that the pairing field A r is smooth on the length 
scale of the relevant fermion frequencies is too crude. As in the case of the Thomas-Fermi 
approximation, which breaks down in the region where the order parameter becomes zero, 
we expect that the assumption that A r varies smoothly breaks down in a region near the 
point where A r becomes zero. In figure 1, the results deemed unreliable are shown in thin 
curves, whereas the results in the region where the aforementioned assumption is estimated 
to hold are shown in thick curves. 
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FIG. 2: The vortex core size, in units 1/&F, is shown as a function of the interaction strength 
l/(&Fa s ). In the BEC regime, the curve shows the molecular healing length, whereas in the BCS 
limit the result goes to a constant as indicated by the line. 



Thus, we must conclude that this assumption may bring us qualitative insight into the 
core region, but that the quantitative analysis of how exactly A r and p r depend on r near 
the vortex core is misleading. It is not yet clear whether the current result obtained by the 
saddle-point approximation can be retrieved in the framework of a Bogoliubov-de Gennes 
theory where similar assumptions (i.e. that A r varies smoothly on the length scale corre- 
sponding to the relevant fermion frequencies) are made. 

Still, the present results allow us to extract a value for the core size, since for this purpose 
one does not need to know the exact functional dependence of A r on r. We estimate the size 
of the vortex core as the distance from the core at which our solution for A r becomes zero. 
The result is shown in figure 2. In the molecular BEC limit, one expects that the core size is 
related to the healing length for the molecular Bose-Einstein condensate, £ m = 1/ \/8im m a m . 
n the path-integral treatment to lowest order in the fluctuations around the saddle point 
^], the molecular scattering length is twice the atomic scattering length, a m = 2a s . The 
density of molecules is half the density of atoms if all atoms form molecules, n m = k F /(6ir 2 ). 
Thus £ m = a/37t/ (8kpa s ) if we express £ m in units of hp 1 . Surprisingly, in the BEC limit 
we find that the vortex core size derived with the path-integral method follows exactly the 
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molecular length £ m . Approaching the crossover region, the vortex core size deviates from 
its molecular value. In the BCS limit, the vortex core size tends to a constant, given by 
l/(2kp), as indicated in figure 2. This suggests that for vortices in ultracold dilute Fermi 
gases the vortex cores in the BCS limit can be much smaller than the BCS correlation length 
£bcs (since &f£bcs ~ E F /A 3> 1), as is the case superfluid vortices in neutron stars 



V. CONCLUSIONS 



The path-integral formalism allows to investigate the ground state of the fermionic su- 
perfluid both in the BCS, BEC and crossover regimes in a unified manner, retrieving the 

nn 

correct limiting behavior [7|, |8||. In this paper we have extended the path- integral treat- 
ment to investigate the vortex state of a fermionic superfluid, by introducing an appropriate 
saddle-point for the pairing field. To set up and solve the resulting saddle point equations, 
the assumption was made that the pairing field varies slowly on the scale of the relevant 
fermion frequencies. Limitations of this assumption become apparent as one tries to calcu- 
late the exact density profile across the vortex core. Nevertheless, with this assumption the 
vortex core size can be calculated as a function of l/(kpa s ). The vortex core size varies in 
the BEC limit according to the molecular BEC healing length, supporting the formalism. 
The path- integral formalism furthermore predicts that the vortex core size (in units of 1/ hp) 
tends to a constant value as l/(kpa s ) is tuned into the BCS regime. 
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